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Abstract 

The vortex-like solutions are studied in the framework of the gauge model of discli- 

nations in elastic continuum. A complete set of model equations with disclination driven 

dislocations taken into account is considered. Within the linear approximation an exact 

solution for a low-angle wedge disclination is found to be independent from the coupling 

constants of the theory. As a result, no additional dimensional characteristics (like the 

CO ' core radius of the defect) are involved. The situation changes drastically for 2tt vortices, 

where two characteristic lengths, and lyy, become of importance. The asymptotical 

| behaviour of the solutions for both singular and nonsingular 2ir vortices is studied. Forces 

O ' between pairs of vortices are calculated. 

ON ' 
ON 

o 

1 Introduction 

Topologically nontrivial objects arising in various physically interesting systems are the subject 
of considerable current interest. It will suffice to mention the 't Hooft-Polyakov monopole in 
^ ■ the non-Abelian Higgs model, instantons in quantum chromodynamics, solitons in the Skyrme 
model, Nielsen-Olesen magnetic vortices in the Abelian Higgs model, etc. (see, e.g., the 
book|IJ). Notice that similar objects are known in condensed matter physics as well. For 
instance, vortices in liquids and liquid crystals, solitons in low-dimensional systems (e.g., in 
magnetics, linear polymers, and organic molecules) as well as the famous Abrikosov magnetic 
vortices in superconductors are the matter of common knowledge. Mathematically, all these 
objects appear in the framework of nonlinear models as partial solutions of strongly nonlinear 
equations. As is well-known, there are still no general methods to study such equations. This 
makes the derivation of the solutions difficult. An important point is also that all the solutions 
are topologically stable and belong to nontrivial homotopic sectors. 

It should be noted that elastic media also leave room for topological defects known as dis- 
locations and disclinations. However, these defects are usually treated within the linear theory 
of elasticity. In this case, all the information about defects is incorporated via source terms in 
the equilibrium condition. The sources are assumed to have a 5-function form multiplied by 
either the Burgers or the Frank vector, respectively (see, e.g., @). There were some attempts 
to invoke the nonlinear theory of elasticity for description of dislocations (see details in ||) 
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which account for the nonlinear relation between stresses and strains thus giving a possibility 
to determine stress and strain fields near a dislocation at large deformations. Recently Q, 
nonlinear problems in dislocation theory were studied in the framework of the gauge model of 
dislocations and disclinations proposed in |5|, |J. Unfortunately, all these approaches || f| are 
based on the perturbation scheme. It is known, however, that topologically nontrivial solutions 
cannot be found in the framework of any perturbation scheme. 

As has been shown in J7], [| |9], [L0|, the gauge model admits exact vortex- like solutions 
for wedge disclinations. This finding confirms the view of a disclination as a vortex of elastic 
medium. Notice that though the correspondence between dislocations and vortices has been 



known for many years (see, e.g., review [J I]) the explicit vortex-like solutions for disclinations 
were obtained only within the gauge-theory approach. It is interesting to note that the elastic 
flux due to rotational defects was found to be completely determined by the gauge vector fields 
associated with disclinations. In continuum there are no restrictions on the value of the flux 
which is in fact the Frank vector, f2. We will consider in the present paper two cases: low-angle 
vortices with |f2| <C 1, and 2tt vortices with |fi| = 2n. 

Before proceeding, let us mention that disclinations are of importance in various crystalline 



and non-crystalline materials (see, e.g., [12]). Among the current applications one can name 



the new class of carbon materials: fullerenes and nanotubes (see, e.g., fll3"| ). According to 
the Euler's theorem, these microcrystals can only be formed by having a total disclination of 
47r. According to the geometry of the hexagonal network it means the presence of twelve 60° 
disclinations on the closed hexatic elastic surface. Notice that disclinations in liquid crystals 
is one of the best-studied cases. In particular, the known exact solution, hedgehog, has been 
obtained within the continuum model of nematics ||. 

It is interesting that a hedgehog-like solution was also found for a point An disclination in 
the framework of the gauge model [[TJ], |T5|. An important advantage of the gauge model follows 
from the fact that it is similar to the known field theory models, first of all to the non-Abelian 
and Abelian Higgs models, where topological objects are studied well. Taking into account 
this similarity we have found two exact static solutions for linear disclinations [7|, |9], [H| which 
will be discussed in this paper. In particular, within the linear scheme our model recovers the 
equilibrium conditions of the standard elasticity theory with a disclination-induced source being 
generated by gauge fields. It should be noticed, however, that these solutions were obtained 
within the restricted model where the dislocation-induced contribution to the Lagrangian was 
neglected. As a matter of fact, this contribution always exists (so-called disclination driven 
dislocations ||). The goal of the present paper is to consider the most general model for linear 
disclinations in elastic continuum involving all main contributions. 

The paper is structured as follows. In section 2 the gauge model of disclinations is presented. 
The Lagrangian describes elastic deformations and self-energy of disclinations as well as a con- 
tribution from the rotational dislocations. A complete set of field equations is formulated for 
the static case. In section 3 the gauge model is studied within the linear approximation. An 
exact vortex-like solution for a low-angle straight wedge disclination is found to be independent 
from the coupling constants of the theory. Forces between vortices as well as dipole configura- 
tions are studied. Topologically stable 2tt vortices are considered in section 4 as the solutions 
of strongly nonlinear equations of the theory. First, we consider an exact singular solution and 
show that the dislocation-induced contribution becomes important. Second, we analyze briefly 
a problem of the nonsingular vortex. In particular, appropriate asymptotic solutions are found. 
Section 5 is devoted to concluding comments. 
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2 Gauge model 



In accordance with the basic assumption of the gauge approach f5|, topological defects can be 
introduced into the Lagrangian of elasticity theory through gauge fields. Namely, the defect- 
free Lagrangian of elasticity theory is invariant under the homogeneous action of the space 
group G = SO (3) > T(3). Let us consider this group as the gauge one and assume that the 
inhomogeneous action of G leaves the initial Lagrangian invariant. In this case, we arrive at 
the Yang-Mills-type theory which contains two different gauge fields associated with rotational 
SO (3) and translational T(3) group, respectively. The model is based on the Yang-Mills min- 
imal coupling arguments. While the initial Lagrangian is chosen to be quadratic in distortion 
fields, the presence of gauge fields makes it strongly nonlinearQ That is why the nonlinear 
relation between stresses and strains occurs from the very beginning. 

For disclinations and rotational dislocations, only 5*0(3) group should be taken into account. 
The problem becomes simpler for rectilinear defects, which are of interest here. Indeed, in this 
case the rotational symmetry becomes broken only in the plane normal to the defect line, and, 
correspondingly, the gauge group reduces to SO (2). Notice that such gauge group is typical 
for models containing vortex-like objects. 

The Lagrangian that is invariant under the inhomogeneous action of the SO (2) gauge group 
takes the following form @ 

£ = C x + C,), + Cw-i (1) 



where 



C x = l -p,B\Bl - l -[X{tiEf + 2/,tr£ 2 ] (2) 



describes the elastic properties of the material while 

^ = -\siD\ b k ac k bd Dl d (3) 

and 

C w = - l -s 2 F ab g ac g bd F cd (4) 

describe defect-induced contributions. Here Eab = B l A B % B — 5ab is the strain tensor, D l ab = 
e'jF a i ) x : ', F a b = d a Wb — dbWai S\ and s 2 are the coupling constants, tiE = Eaa and summation 
over repeated indices is assumed. In accordance with the minimal replacement arguments, the 
distortion tensor is written as 

B a = d aX l + e)x j W a (5) 

where x % ( xa ) = X % ( xA iT) characterizes the configuration at time T in terms of the coordinate 
cover (x A ) of a reference configuration, W a is the compensating gauge field associated with the 
disclination field. In (3) the tensor k ab is given by k AB = —5 AB , k 33 — 1/y, and k ab = for 
a 7^ b, whereas in (4) g AB = —5 AB , g 33 = l/£, and g ab = for a ^ b. The parameters y and £ 
are the two positive "propagation parameters", e* is a completely antisymmetric tensor, e\ — 1, 
and A and fi are the Lame constants. (Indices a,b,c,... =1,2,3 and A,B,C,...=1,2 are the space 
labels, whereas i,j,k,...=l,2 belong to SO (2)). 

Notice that Cw describes a self-energy of pure rotational defects (disclinations). In accor- 
dance with the general approach it acquires a standard form of the SO(2) Yang-Mills action. L$ 

1 The cubic terms can be included into the Lagrangian as well (see 0), but in so doing the model becomes 
too cumbersome for analysis. 
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is an additional invariant term associated with translational defects (rotational dislocations). 
This is reflected in the fact that the tensor D l ab in (3) is directly related to the density of 
dislocations defined as a Al = e ABC D BC (see ||). Thus the appearance of this term in the La- 
grangian reflects the known in the defect theory fact that the presence of disclinations implies 
the presence of dislocations. 

The Euler-Lagrange equations for (1) take the following form in the static case 

d A a A = e{W c af + s lX l 8 lk F AB F AB , (6) 

a A 6 l3 elx l = 2d B {(s 1 ( X ) 2 + s 2 )F BA }, (7) 

where (x) 2 = X l Xu Z = 1,2. To avoid cumbersome expressions we will sometimes omit the 
right order of the top and bottom indices which can be easily restored by using the appropriate 
(5-symbols. The stress tensor cr J c is determined to be 

a^ c = i[A(trJ5)fl£ + 2^E CB B j B )]. (8) 



It is convenient to introduce the dimensionless variables via x % = Js 2 /s\X l and W A 



Si/s 2 Wa- The Euler-Lagrange equations (6) and (7) become 



d A a A = e{W c af + ^x^ikF AB F ABl (9) 

&A^Ai? = —Um 2 + i)F B % (io) 

HS 2 



where x\ xB ) — \J s 2/ Six\ xB )i d B = \Jsi/s 2 d B , F AB = d A W B — d B W A and the stress tensor is 
found to be 

af = ^(trEI)& A + E A c& c . (11) 
Here L = X/fi, and the strain tensor takes the form 

E AB = B A B B — Sab (12) 

with B\ = d A x l + CjX'WA- To simplify notation, below we will omit the "tilde" symbol. As is 
seen, the self-consistent system of field equations (9) and (10) is strongly nonlinear. In some 
respects (see discussion in 0, [|) it is similar to that in the Abelian Higgs model. An analogous 
system of equations appears under study of type-II superconductors in magnetic field directed 
along the z-axis. This observation was helpful in finding of the exact vortex-like solution to 
(9) and (10) for si = 0. A possible way to study these equations is the linearization 
procedure which is valid for low-angle defects. It was found |10| that there is an exact vortex- 



like solution at si = which reproduces the known strain and stress fields of straight wedge 
disclinations. Notice that both solutions are found to be singular at the disclination line. While 
this singularity is well known in dislocation theory, the question arises if there exist nonsingular 
solutions. Besides, it turns out that for Sx = the coupling constant s 2 drops out of the problem 
as well. Let us reexamine the vortex-like solutions within the most general model. 
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3 Low-angle disclination vortices 



Let us consider first the linearized equations of the model. The linearization procedure is based 
on a homogeneous scaling of the gauge group generators (see details in |5|]). It is clear that 
employing the linear approximation we restrict our consideration to small deformations or, 
correspondingly, to the low-angle defects. The displacement vector u l can be introduced as 
follows 

X \x B ) = 5\x A + u\x B ). (13) 
Then, with the scaling parameter e all the fields are expanded in series of e 

v* = eu{ + e 2 u{ + W A = eW 1A + e 2 W 2A + ... (14) 

Taking into account only defect-induced displacements u\ that are of interest here, we get the 
first order equations in the following form 

d B [(x 2 + y 2 + l)F BA } = 0, (15) 

s 2 

Au + (L + l)Vdivu = j + x—F AB F AB (16) 

where the density of elastic flow due to a disclination, j, is completely determined by gauge 
fields 

jc = -{L - l)t AC W A - Le AB x B d c W A - e AB x B d A W c - t C BX B d A W A . (17) 

Hereafter we omit the index 1 denoting the order of the approximation. Let us emphasize that 
we assume here sf///s 2 ~ l/ e - F° r other two possibilities, s\/ /is 2 ~ e and s 2 /fis 2 ~ 1, terms 
with Si in (15) and (16) are of little importance thus the standard theory |TU| is recovered. 
Notice that for si = a solution of these equations was found in flQ|| . An interesting property 
of this solution is its independence from the parameter s 2 as well. This can be seen directly 
from (15) and (16) where s 2 is completely absent. Let us choose the following vortex-like ansatz 

W A = -W(r)e AB d B logr, u* = x l G{r) (18) 

where r 2 = x 2 + y 2 . With (18) taken into account one can rewrite (15) and (16) as follows 

W'(r 



r 



r 2 + l)- 



(19) 



(L + 2) ^(r) + ^= L^-^ + M (m) a (20) 
\ r J r r 2 fxs 2 \ r J 

where G' stands for dG/dr and W for dW/dr. A solution to (19) takes the form 

W(r) = Ct ln(r 2 + 1) + C (21) 

with Cq and C\ being the integration constants. Notice that a disclination flow through the 
plane xy is given by a circular integral 

^- j> Wdr = v. (22) 
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Taking this into account, we immediately obtain that G\ = 0. Thus the constant Cq turns out 
to be in fact a topological characteristic of the defect, that is the Frank index v. For W(r) — v 
(20) becomes remarkably simpler and has a solution 

G(r) = -j^ lnr - + C 3 . (23) 

Since the boundary condition requires tt l (0) = we must put C 2 = 0. Turning back to the 
dimensional variables, we finally obtain 



u i = -x i (-—lnJ^r + C 3 ), (24) 




where C3 is still an arbitrary constant. As is seen, the term with s x and s 2 only renormalizes 
the constant C3. As an example, for the straight wedge disclination on a disk of radius R with 
a boundary condition in the form u l {R) = one obtains 

tf = -'rT2 ln s- (25) 

Similarly, for the most-used boundary condition, <Jki n i = at the free surfaces, one can repro- 
duce the well known stress fields for a wedge disclination on a disk (see details in [[Hj). It is seen 
that parameters s\ and s 2 actually drop out from (24) and (25). Thus one can conclude that the 
information carried by the coupling constants s\ and s 2 is lost within the linear approximation. 
What does it mean? As is known, the classical theory of elasticity introduces a characteristic 



velocity Jfi/po, but does not lead to a characteristic length. For this reason there is no room 
within the linear theory of elasticity for description of the core region. It is interesting that the 
gauge theory of defects || |5| introduces appropriate length scales. These are the dislocation 
length scale, ti = s\f [A, and the disclination length scale, % = s 2 //z. Nevertheless, as we have 
just seen, in the linear approximation the gauge theory loses these parameters thus making 
the description of the core region impossible. One can expect^, however, that these parameters 
would be of importance in a study of the basic model equations (6) and (7). We will consider 
this problem in section 4. 



3.1 Forces between vortices, dipole configurations 

Let us consider two low-angle vortices with parallel Frank vectors oriented along the z-axis. This 
corresponds to a pair of straight wedge disclinations. For simplicity, we suppose that disclination 
lines coincide with their axes of rotations. The stress field due to the disclination results in 
the force acting on the second defect (in perfect analogy to the known Peach-Koehler force in 
dislocation theory). Generally it can be written as [M (per unit length of the disclination line) 

•F z ^bkc^amn^m-^-n^ab^ki (^6) 

— * 

where £ is a tangent vector at the disclination line, Vt m are the components of the Frank vector, 
X n = x n — and x° is a point on the axis of the disclination. In our case, = and one 
has to put in (26) k = 3 and m = 3. Notice that the same expression follows from the general 
equations of the gauge model |5|, |J in the linear approximation. As was shown above, the stress 

2 We would like to thank Prof. A.M. Kosevich for attracting our attention to this possibility. 
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fields take the well-known in the linear theory form. Evidently, the force between two parallel 
low-angle wedge disclinations also has the known form. For example, when the first vortex is 
situated at point (0,0) while the second one at point (d,0) on the plane xy we obtain from (26) 



T x = dQcr. 



yy ■ 



T 



-dQa 



tj.i; 



(27) 



in accordance with [[Uj. Notice that for vortices with equal but oppositely directed Frank 
vectors such configuration corresponds to a wedge disclination dipole with nonskew axes of 
rotations. It is interesting to reproduce the solution for the dipole within the gauge model. 
Since the previous analysis shows that the constants s± and s 2 are inessential in the linear 
approximation, we will drop terms with s± and put W(r) = v from the beginning. A dipole 
solution to (15) then reads 



W B = -ve B cdc(logn - logr 2 ). 



(28) 



The most simple way to solve (16) is via the Airy stress function. Namely, let us differentiate 
(16). After straightforward calculations one can rewrite this equation as 



(A/4/x + 1/2) AtrE = e AB d A W B . 



(29) 



The last term in the right-hand side of describes a source due to disclination fields. For 



solution ( pq) it takes the form 

e AB d A W B = ^A(logrx - logr 2 ) = 2itv{8(r x ) - 6(f 2 )). 

Introducing the Airy stress function x(r) by <jba = ^bm^an9m9nx{^) ■> one can finally rewrite 
the equation (|29|) as 

K^ 1 A 2 x = 2^(5(f 1 )-5(f 2 )). (30) 

Here K = 4/i(A + /i)/(A + 2/i), and txE = (1/(A + fi))Ax(r). Evidently, a solution to (|30|) 
is the sum \ = xi + X2 with x% = Arflnri (i=l,2). One can easily find that A = dzz/i^o/4. 
Finally, turning back to oba one can exactly reproduce the known stress fields for disclination 



dipoles (cf., e.g., |L2], |T 



4 2tt disclination vortices 

Let us choose the following ansatz for (6) and (7) to meet the necessary symmetry requirements 

X X {x A ) = F(r)cos#, x 2 (x A ) = F(r)sm6 (31) 

and 



WJx J 



y 



W(r) 



Wy{X t 



X 



W{r) 



(32) 



where r 2 = x A x A (r, 6 are the polar coordinates). All the variables here are again dimensionless. 
We restrict ourselves by the topological sector with n = 1. As is seen, (31) and (32) describe a 
27r vortex, that is the circular integral in (22) is equal to v — 1. With (31) and (32) taken into 
account (9) and (10) reduce to 



F 
4^W 



W 



f K{g 2 + f - 2) + 2P{f - 1) + K 



T (fg) + — 

dr r 
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+ 



3(K + 2P)g'g 2 
. d 



dr 



1 + F' 



-2{K 
W'~ 



P)g' 



1 r 



+ - 

r 



(K + 2P)g 3 -2(K + P)g 



Ff[K(f 2 + g 2 -2) + 2P(f 2 -l)] 



(33) 
(34) 



where K = \s 2 /sj, P = [is 2 /s 2 , g = dF(r)/dr, f = F(r)(l - W(r))/r, W = dW{r)/dr. This 
system of equations is of our interest in this section. We will consider two possible cases. 



4.1 Singular vortex 

In the case of s\ = an exact solution to (33) and (34) for a static disclination vortex was 
found in ||. The solution is singular on the disclination line with W(r) = 1. It is interesting 
to note that the same solution is valid for the general case when s\ ^ 0. Indeed, for W(r) = 1 
one gets / = 0. In such an event, both sides of (34) turn out to be zero whereas (33) reduces 
to 

(3g'g 2 - N 2 g') + i (V - N 2 g) = (35) 

where iV 2 = 2(K + P)/{K + 2P) = 2(A + /i)/(A + 2fi). Carrying out an integration in (35) one 
obtains finally the algebraic equation 

\g*(r)-N 2 g(r)\ = ^ (36) 

where Cq is an arbitrary integration constant. The solution to (36) is written as g(r) = 
(2/VS)N g(r) with 

-/ ^_ / 9i( r ) = cosh[|cosh -1 (r /r)], r < r 

g{r) — < . , [61) 

{g2{r) = -cos[g cos L {r /r) + — J, r > r 

where r = 3y/3C /2N$ is a characteristic parameter, and / = 0, 1, 2. Notice that in the natural 
variables ro becomes dimensional. It was supposed in that ro could be considered as a core 
radius of the defect. Indeed, according to (37) the point r = ro is prominent. It should be 
mentioned, however, that this attractive possibility for description of the core radius requires 
an involvement of the additional phenomenological parameter, Co, into the theory. At the same 
time, the model parameters Sj and s 2 drop out of (35) and (37). 

The main reason is that the chosen ansatz W(r) = 1 (the pure gauge for all r) is too 
restrictive. Obviously, any solutions with no constant W(r) are of special interest. However, 
equations (33) and (34) look rather cumbersome, and a search for nontrivial exact solutions 
remains still an open problem. Let us try another way of looking at the problem. For this 
purpose, one can put C = in (37). The exact solution to (35) takes then the essentially 
simpler form 

W(r) = l, F lj2 (r) = ±iV r, F 3 (r) = 0. (38) 

As a first step, let us consider small perturbations of the exact solution (38). It will be shown 
below that even this simplified consideration allows us to get an important information about 
the role of coupling constants si and s 2 in the theory. Let us write 

F(r) = N r + eui(r) + e 2 u 2 (r) + (39) 

W(r) = 1 - ewi(r) + (40) 
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We consider here the case F(r) = Fi(r). Substituting this expansion into (33) and (34) one 
gets 

(i + iV 2 r) w'[{r) + (iV 2 - 1) w[(r) - ^N*r Wl (r) = 0, (41) 

It should be noted that U\ = 0. This follows from the requirement F(r) — > at r — > 0. Let us 
analyze (41). As is known, the equation of the type w" + H(x)w' + Q(x)w = can be put into 
the form z" + I(x)z = by a substitution 



w(x) = z(x) exp ^— - J H{x)dx S j 



where J = -\H' - \H 2 + Q. For (41) 



N 2 r 2 - 1 „ _ PN, 



i r 2 



*< r > = nk^ry Q{r) = - wfey (43) 



In this case, the general solution is not known yet. Instead, let us derive two limiting cases. 

Lin the limit N^r 2 3> 1 one obtains H(r) = 1/r, and Wi(r) = z(r)/y/r. The equation for 
z(r) takes the form 

'PN 2 1 + 2P N 
~2 4^ 



Ar)=(^~^U(r). (44) 



This is a special form of the Whittaker equation with the solution 

z(r) = QW^ m {f3r) (45) 

where C\ is a constant, Wo^ m is the Whittaker function, m = ±i^J P/2, and fi = N \ / 2P. In 
this case, W(r) is found to be 

W(r) = l-C l r- 1 / 2 W , m (M- (46) 

Notice that C\ includes e. Depending on fi two asymptotics can be obtained. 

(i) For fir 3> 1 one gets 

W(r) = 1 - C,r _1/2 exp(-i/3r) (47) 

and 

F M = ^ + c 'i(ra)A^ ex P(-«- < 48 > 

(ii) In the limit /3r <C 1 one has 

W(r) = 1 - d cos (/p/2m(/3r)) (49) 



and 



F(r) = ^(2 - Cf)r. (50) 
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2. Let us consider the limit N^r 2 < 1. In this case, H(r) = —1/r and W\(r) = \/rz(r). 
The equation for z(r) reads 

z"(r) = (A + ^iv V) *( r ) (51) 



with a solution in the form 



, , sinh(l7 2 r 2 ) , 
z(r) = 2C S (52) 



where 7 2 = N$JP/8. Thus 



r 



2(7 1 

W(r) = 1 - — sinh(- 7 V)- (53) 



77 v 2 



Two limiting cases are of interest, 
(i) For 7 2 r 2 3> 1 one obtains 



v/7 v 2 

and 



C 1 

W(r) = 1 - exp(- 7 V) (54) 



F W = iVor + ^3|(7^_l_ (M) 



where ^ = iV^v 7 ^ - AT)/(AT + 2P). 
(ii) If 7 2 r 2 <lwe get 



and 



W(r) = 1 - C s Y /2 r 2 (56) 

F(r) = N r + ^V^^ 3 . (57) 

For better understanding of the obtained results let us return to the dimensional variables. It 
is interesting that both characteristic lengths of the theory turn out to be involved. Namely, 



the dimensional r reads r = y^/sir, and the important parameters (3 and 7 are determined 

via y/P = Iw/ll- m three physically interesting limits were discussed. In particular, in 
accordance with the typical condition which is valid in crystals and polycrystals is lw ^> l<f>, 
in some polycrystals and amorphous bodies lw ~ I4, while the most exotic limit which can 
be expected in some special amorphous materials is lw "C /</,. Our consideration shows that 
vortex- like solutions have different asymptotics in each case. For lw 3> l</>, lw ~ l<t>i an d lw ^ l<j> 
they follow (1 (i) ,2 (i) ) , (l(i),2(ii)), and ( 1 (ii) ,2 (ii) ) , respectively (see above). In accordance with 
(11) and (12) this results in different asymptotic behaviour of strains and stresses due to 
disclinations thus giving a possibility for the experimental verification. On the other hand, 
the obtained results indicate that the proper information about the core region, if any, can 
be obtained only within the framework of the complete gauge model which should include 
rotational dislocations. 



4.2 Nonsingular 2tt vortex 

Let us discuss briefly a possibility for a nonsingular solution in (31) and (32) which provides a 
finite energy of the vortex. This means that condition W(r) — ► for r — > should be satisfied. 
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A simple analysis of (33) and (34) shows that there is an asymptotic solution at small r in the 
form 



W(r) ~ r a , 



F(r) ~ ar A 



(5* 



where /i = 1, a = 2, and a is an arbitrary constant with the only restriction a ^ 1 following 
from (34). Notice that this resembles the behaviour of the Abrikosov-Nielsen-Olesen vortex. 
For large r, the asymptotics found in the previous subsection are valid. To prove the existence 
of the solution for any r the numerical calculations of variations in the energy density which 
has the following form 



E(r) 



A 



+ - 

S2 



f + — (i - wy - 2 



+ 2/x 



f + — (i - wy - if' 1 - 2— (i - wy + 2 



(59) 



where i 7 " = dF/dr and W = dW/dr must be used. This study will give a final answer about 
whether the above asymptotics come from the unique solution or not. The corresponding 
calculations are now in progress. 



4.3 Forces 

Let us discuss briefly the force between two 27r vortices. In accordance with the classical formula 
(26) one obtains 



T A = Qf X(g 2 + f- 2) + 2^-1] 



x A . 



(60) 



Thus the force turns out to be exactly zero for the solution W(r) = l,F(r) = N r, that is in 
the case of 2tt singular vortices. Assuming a small perturbation of this exact solution we get 
the following asymptotics: 
1. For N$r 2 s 1 /s 2 > 1 



-1/2 



T A = -2^Nldn r. 



(61) 



2. A^ 2 r 2 si/s 2 < 1 



T A = -4/iiV 2 C s fi— sinh f 7 2 — r 2 ) x A (62) 
Vl V s 2 J 

It is important to note that in the strict sense we have to use the general expressions for the 
forces given by the gauge model (see @, ||). These calculations, however, are too cumbersome 
and will be omitted here. Notice only that the main conclusions agree very closely with (60)- 
(62). 



5 Conclusion 

In this paper we have studied the vortex-like solutions for disclinations within the most general 
gauge model of rotational defects when both disclinations and rotational dislocations are taken 
into account. The model contains two additional parameters, coupling constants s\ and s 2 , 
which allow us to introduce two characteristic lengths 1$ and lw The appearance of these 
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lengths is the unique property of the gauge model that distinguishes it from the classical 
elasticity theory as well as from other known models of elastic continuum with topological 
defects. 

There are two distinctive features of the vortices in elastic media. First, the elastic flux is 
'classical' in its origin, i.e. there is no quantization as opposed to the magnetic vortex. This 
means that generally there are no restrictions on the value of v in (22) apart from v > — 1 for 
topological reasons. However, if we take into account the symmetry group of the underlying 
crystal lattice the available values of v become 'quantized' in accordance with this group (e.g., 
v = 1/6,1/4,1/3,... for hexagonal lattice). Second, the singular character of the solution 
on the defect line is typical for the dislocation theory. As a result, all the known solutions 
for dislocations contain a logarithmic divergence in the energy. To avoid this difficulty, one 
introduces a cut-off from below by using r as a core radius of the defect. The core region 
itself is assumed to be beyond the scope of the linear theory of elasticity. For this reason, any 
nonsingular solution will be of essential interest. 

Finally, let us note that a similar problem appears for point 4-7T disclinations. In this case, the 
gauge group 5*0(3) should be considered. We expect that the inclusion of rotational dislocations 
will clarify the role of the characteristic lengths as well as the problem of the core region in this 
case. This study is now in progress. 

This work has been supported by the Russian Foundation for Basic Research under grant 
No. 97-02-16623, and the Slovak Scientific Grant Agency, Grant No. 7043. 
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